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Abstract of EP1 093231 

The LDPC (Low Density Parity Check ) coding procedure 
codes M information symbols with a N-M redundant symbols 
with a checking matrix A which has the same minimum 
number of non zero elements in each line and less than two 
lines or columns with only one non zero value. 
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(54) Procede de construction et de codage simple et systematique de codes Ldpc 



(57) Le procede permet la construction de codes 
LDPC comportant N symboles dont K libres chaque co- 
de e tant defini par une mathce de contrdle A comportant 
M=N-K Hgnes N colonnes et t symboles non nuls dans 
chaque colonne. II consiste: 

- & attribuer a toutes les lignes de la matrice de con- 
trdle A un merne nombre de symboles non nuls. 

- & prendre comme nombre de symboles t un nombre 



impair le plus petit possible. 
k definir les colonnes de facon que deux colonnes 
quelconques de la matrice de controle A n'aient au 
plus qu'une valeur commune non nulle, 
et a definir les lignes de facon que deux lignes de 
Ja matrice de contrdle A n'aient qu'une valeur com- 
mune non nulle. 

Application: Transmissions radio. 
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[0020] Le procecte de construction est decrit ci - apres selon deux variantes, selon que rn est ega! a 1 ou t. 
[0021] Les valeurs difterentes de m ne sont pas a envisager ici. a cause de la condition "a u qui exige que tr soit 
entier. En effet, tr = t N / M, soit encore t n / m. 

[0022] n et m etant premiers entre eux, m doit diviser t et ne peut done £tre qu'Sgal a 1 ou t pour t premier et petit 

5 (vrai pour les faibles valeurs de t, a savoir 3, 5, 7). 

[0023] Dans la premiere variante oum=i (codes de redondance r/ (r - i))je procede selon Invention est valable 
pour les codes de taux de redondance N/K de la forme r / (r - 1), ou le nornbre de syrnboles de redondance est 
exactement N / r (r entier). Dans ce cas, M est egal a P et le tableau de la figure 2 se resume au tableau de la figure 
3. Le precede" consiste alors a chercher n sequences de longueur M comportant t "1" et (M - 1) "0". 

10 [0024] Ces sequences not^es ci-apres w[0..n-i], sont obtenues par : 

- une autocorrelation cyclique 6gale a 0, 1 , ou t (la sequence i decalee ne coincide avec elle-meme non d£cal<§e 
qu'en 0 ou 1 point) telle que pour tout i = o...n-l : 

15 

z jm>..../w-i wTOM wfflk] = t (par definition) 



£k=o..a4-i w[i][k]w[i][(k+p) modulo M] = 0 ou 1, pourp= 1..M-1 

20 

- une intercorreiation cyclique egale a 0 ou 1 (la sequence i decalee ou non ne coincide avec la sequence j qu'en 
0 ou 1 point) telle que pour toute paire {i. j} ou i = 0...n-l et j - 0...n-l sont differents : 

25 i k-o....am w BJM w 0l[(k + P) modulo M] = 0 ou 1 , pour p = 0..M-1 

[0025] L'algorithme de calcul des sequences w est tres simple : il determine successivement ies positions pos[0][0... 

t-1], pos[l][0...M] pos[n-i][o...MJ. ou ces sequences possedent un T, en commen-ant parpos[xp] = 0, pos[x] 

f 1 J = 1 pos[x][M ) = t-1 , et en les modifiant pour satisfaire les conditions d'autocorr6lation et interco relation. 

30 [0026] Pour t a 3, les boucles de calcul mises en oeuvre sont montr^es a I'annexe 1 . 

[0027] Cet algorithme 6choue lorsque n est trop grand compte tenu de M : peu de "petits" codes sont trouves. mais 
ceci est de peu d'importance car on cherche gendralement des codes de grande taille (N » 100). 
[0028] Les colonnes de la matrice A sont alors tout simplement les vecteurs w permutes circulairement : 

35 - kerne sous-matrice (k = 0...n-1 ) 

A[ligne][colonne) = w[k][(ligne - (colonne - k P)) modulo M] 
Avec: 

ligne = 0...M-1 

et colonne = kP...(k+ 1)P - 1 

40 

[0029] Ainsi. chaque ligne de A comporte exactement t valeurs non nulles dans chacune des n sous-matrices soit 
un total de n = tr. 

[0030] Un exemple de matrice A. obtenue par ce proc6d§ pour un code LDPC(75, 50. t=3, tr=9) de redondance 3 / 
2 (r = n = 3) avec P = 25 est montre a la figure 4. Selon le tableau represents on peut constater que: 

45 

w[0]p]= 1 pouri= 0, 1, 3 
w[1][ij= 1 pouri = 0, 4. 9 
w[2][i] = 1 pour = 0, 6,13 

so [0031] La construction proposee garantit que : 

- chaque colonne comporte exactement t valeurs non nulles (par definition des w) 

- chaque ligne comporte exactement tr valeurs non nulles (grace aux proprtetes ^autocorrelation et intercorrSlation 
des w) 

toute paire de colonnes distinctes a au maximum une valeur non nulle commune (idem) 
toute paire de lignes distinctes a au maximum une valeur non nulle commune (idem) 

[0032] Selon une deuxieme variante inspiree de la precSdente correspondant au cas ou m=t le procede selon I'in- 
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c ff-o....Afi w t°][ k ] w[0][k] - t (par definition) 



z k=o....am w [0J[k] wOJ[(k + p) modulo M] = 0 otM , pour p =- 1 ..M-1 



[0034] En fait, c'est la m§me definition que pour m = 1 
[0035] Les sequences suivantes, w[i ...n-m] sont obtenues par 

pour tout i= i...n-m: 

z ff=o ...M-i w PlW w[i][k] = t (par definition) 
£ k=o.../w-i wf'JM w[ij[(k + p m) modulo M] = 0, pour p = 1 ..P-1 

pour tout i= 1...n-m : 

z k=o...m-i HOW w[0][(k + p) modulo M] = o ou 1 , pour p = 0..M-1 
• et par une interaction cyclique avec les sequences w[l..n-m] egale k 0 ou 1 pour des decalaaes m„ir.m« h* 

l /c=o....m.i w PJ[k] wtj][(k fpm) modulo M] = 0 ou 1 , pour p - 0..P-1 

Cq ™ 68 " 8 * maWe " S< "" 810,8 " s " M " pe, ™ ,6s ***•««* a«« m pas Pgal a , ou m 

^5 - Sous-malrice M x M k gauche de A : 

A[ligne][colonne] = w[0][(ligne - colonne) modulo Ml 
Avec: 



50 



Hgne = O...M-1 
colonne = 0...M - 1 



- Sous-matrices MxP suivantes (en nombre egal an-m) pour k = m n-1 
A[l.gne][colonne] = w[k - m +- l][(ligne - m (colonne - k P)) modulo Mj 
Avec. 

ligne = o...M-i 
55 colonne = kP...(k + l)P-i 
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[0039] Un exemple de matrice A, obtenu selon la deuxieme variante du procede selon invention est montre & la 
figure 5 pour un code LDPC(75. 30. t=3. tr=5) de redondance 5 12 (n = 5. rn = 3) avec P = 15. On constate sur le 
tableau de ia figure 5 que: 

5 wfojfi] = 1 pour = 0, 1,3 

w[1][i]= 1 pouri = 0,4, 8 
w[2][i]= 1 pour = 0,5, 10 

[0040] Le procede selon ^invention qui vient d'etre decritsous ses deux variantes conduit directement a un algorithme 
10 de codage de symboles de redondance Yi et de symboles d'information Xi tres simple. 

[0041] Pour cela, il suffit de considerer que les symboles de redondance Yi sont les M premiers symboles d'un mot 
de code, et les symboles iibres Xi (information) sont les N - M derniers. 

[0042] Les equations que doivent verifier tout mot de code peuvent done §tre reecrites sous la forme : 

15 

^/-m...a/-i Am ' Xi = 0. pour m=0...M-l 

ou encore : 

20 

z /=o...m-i Ar ™ Yi = Zm , pour m=0...M-1 

avec 



Zm = - ^ hM NA Ami Xi, pour m=0...M-1 

[0043] Le procede consiste alors & calculer dans un premier temps les M quantites Zm de la matrice de passage, 
et ensuite les symboles de redondance : 

30 

Ym = l /=0 /WL1 Bmi Zi, pour m=0...M-1 

[0044] A litre d'exemple, pour le code LDPC(75, 50) les quantites Zm sont calcules par le systeme d'6quation defini 
par le tableau de la figure 6 qui, apres resolution . se transforme en le tableau des symboles de redondance de la figure 7. 
[0045] La matrice B d'element generique Bij est I'inverse de la partie gauche A^ (de dimension M x M) de la matrice 
A. Eile a une forme tres simple : par construction, toutes ses colonnes sont des permutations circulaires de la sequence 
w[0][0..M-1J: 

Aij = w[0][(i - j) modulo M], i=0...M-1, j=0...M-1 
[0046] La matrice B comporte alors M lignes qui sont des permutations circulaires d'une ligne unique b[0...M-l], a 
savoir: 

Bij = b[(j - i) modulo M] 

B etant I'inverse de A M , les coefficients b sont definis par : 

z m>...am w [°]['l b[(i + k) modulo M] = 1 si k = 0. 0 si k = 1 ...M-1 

[0047] Par exemple, pour le code LDPC(75, 50), les coefficients de redondance Ym sont calcules par le systeme 
d'equations defini par le tableau de la figure 8 qui apres resolution se transforme en le tableau de la figure 9. 
so [0048] Cependant, il existe des cas ou le calcul est impossible. 

[0049] On peut en effot ecrire les equations qui les dSfinissont sous la forme : 

*A M t{b[0], b[l] b[M-1)} = t{1,0, 0 0} 

[0050] La matrice ^ est une matrice circulante, dont ia premiere ligne est egale a a[O...M-i] = w[0]. 
[0051] Son determinant est egal au produit de ses M vateurs propres I0...M-1 . 
55 [0052] La kerne valeur propre est elle-meme donnee par : 



5DOCID: <EP 1093231A1 J_> 



5 



EP 1 093 231 A1 



ou « est une racine Meme de I'unitt 
[0053] Parexemple: 



pour w[0] = a[O...M-1J= {1,1,0,1 0 0 0 } 



on a : 



15 



*i =1 + a + a 3 
[0054] Si M est multiple de 7, il se trouve que I'equation 



1 + a + a 3 = 0 



I'on garde ce w[0] ci car : *" ' a,9 ° nthme pn ^' » <™ clone eliminer ceux ou M est multiple de 7 si 



55 



rune des valeurs proprcs de 'A M sera nulle 

done son determinant sera nul 

done or, ne pourra pas trouver de bfi] convenables 

done on ne pourra pas effectuer le codage (caiculer fcs Yi) 



40 



3g Sne"^ 

[0057] On montre facilement (en factorisant x« - 1 e t afx> - v an k, 

d une valeur MO pour laquelle a(x) divise xMO - 1 ~ 1 ~ m W X ) que Ces valeurs de M s ° n < les multiples 

[0058] Par exemple. pour des codes binaires avec t = 3: 

• w[0]={1, 1,0. 1....} 
35 - w[0]={1,o, 1, 1, ...} 

interdisen. M multiple de 7 (a(x, definit une racine 7#me de runite) 

• w[0]= {1,1, 0,0.1,...} 
- w[0] = (1,0, 0,1,1,...} 

interdisent M multiple de 15 (a(x) definit une racine I5eme de funite) 

: m = li: ?: 1 o°: i ;!T pas accept * (autocorr6,ation 

" w[0]= (1,0,0, 0, 1, 1, ...} 

interdisent M multiple de 3 fa(x) est mnitinio ho 1 

(a ( x) est multiple de 1 + x * x2 qui definit une racine 3eme de ('unite) 

" wtOJ={1,0, 1.0,0, 1,...} 
" w[0]= {1,0,0, 1,0, 1,...] 
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est formee de la concatenation de sequences maximales (Maximal Length Sequences) de longueur 7 (ou 15 ou 31 
ou 63). 

[0061] On calcule done la reponse impulsionnelle du filtre a reponse impulsionnelle infinie (Rll) 1 / A(z) et on en 
extrait une tranche de longueur M qui. une fois periodisee, donne la s§rie {1 . 0. 0. ...} apres filtrage par A(z). 
5 [0062] Par exemple. pour un code binaire avec t = 3 et pour lequel seuls a[0], a[kl| et afk2] ne sont pas nuls, un 
algorithme correspondant est fourni a I'annexe 2. 

[0063] Dans un esprit de simplification pour ne pas effectuer le calcul precedent a chaque codage, ralgorithme de 
codage peut encore etre defini par les k2 derniers elements b[M-K2...M-l] qui, par recurrence (filtrage par 1 / A(z)) 
permettent de recalculer tous les autres. 

10 [0064] Egalement comme la deuxieme phase de ralgorithme de codage standard (calcul des Y a partir des Z) com- 
porte en moyenne M 2 / 2 operations, ce qui peut devenir important pour des codes consequents : la complexity etant 
une fonction quadratique de la taille, et que de plus, il est necessaire de stocker le tableau intermediate 2 (M Elements) 
et connaitre le tableau b (M elements aussi) s'il n'estpas calcuie surplace, cette partie de ralgorithme peut §tre modifiee 
pour n'utiliser que deux tableaux intermediates de trfcs petite taille en reecrivant ainsi les equations donnant les Y de 

is ia facon montree par le tableau de la figure 9 (exemple pour un code LDPC(75, 50)) : 

[0065] Les M - k2 (pour t = 3) premieres lignes sont les M - k2 dernieres lignes du systeme ^equations donnant Y, 
avant resolution. 

[0066] Les k2 demises lignes sont les k2 demises lignes du systeme d'equations donnant Y, apres resolution. 

[0067] II suffil alors de calculer les Y dans I'ordre inverse, a savoir Y[M-1], Y[M-2] Y[0]. 

20 [0068] Le nombre derations a effectuer est alors en moyenne de k2 M / 2 (calcul de Y[M-1] ... Y[M-k2]) suivi de 
t(M - k2) (calcul de tous les autres) soit approximativement (t + k2/2) M : la complexite n'est plus qu'une fonction lineaire 
de la taille. 

[0069] L'algorithme utilise X[M...N] en entree. 

[0070] La partic basso de X (X[0...M-1]) est utilisec comme stockagc temporaire pour les Z : X[0...M-1] stocke Z 
25 [k2...M-i , 0...k2-l] pour eviter un decalage circulaire en phase finale. 

[0071] Les b[i] sont calcuies iterativement surplace, a partir des b[M-k2...M-l]. 
[0072] Le code est defini par deux tableaux : 

- le tableau finB[0...k2-1] des k2 derniers elements de b 

30 - le tableau pos[0...(n - m + 1 )t] contenant les positions des elements non nuls des sequences w[0], w[1J w[n-m]. 

On utilise deux buffers internes de taille k2 : 

- reg[0...k2-l ] pour calculer les b[i] 

35 - temp[0...k2-l] pour stocker les valeurs intermediaires de Y[M-k2...M-1]. 

[0073] L'algorithme complet de codage rapide est alors celui montre a I'annexe 3. 

[0074] Ces algorithmes sont trfes simples a mettre en oeuvre. lis ont entre autre la caracteristique de definir un code 
par tres peu de parametres, a savoir les (n - m + 1)(t - 1) positions non nulles des "1" dans les sequences w et 
40 eventuellement k2 coefficients de codage. M§me s'ils ne donnent pas tous les codes possibles repondant aux condi- 
tions a-d (par exemple, pas lecode (150, 125) de redondance 6/5, qui necessite n=6 sequences wde longueur P=25), 
ils en donnent suffisamment pour que, dans une application quelconque ou N et K sont definis a priori, on puisse 
trouver soit 

45 - un code (NLDPC, KLDPC) avec NLDPC = N et KLDPC = K 

- un code voisin (NLDPC+d. KLDPC+d), avec d faible, qui sera raccourci par non transmission de d symboles utiles 
mis arbiirairement a zero. 

[0075] Par exemple, pourobtenir uncodede (N, K) redondance 5/3 (taux 0.6), ilsuffitde partir d'un code (NLDPC+d, 
so KLDPC+d) de redondance 8 / 5 (taux 0.625) avec d - NLDPC / 1 5. Pour des valeurs de N inferieures ou egales a 500 
et t = 3, il est possible de construire tr£s rapidement 932 codes differents dont les redondances sont les suivantes (on 

s'est volontairement limite aux redondances comprises entre 4 et 8/7, et aux codes ou w[0] = {1, 1, 0, 1 , 0, 0, 0 } 

pour lesquels k2 = 3) : 



R = 


4/1 


soit 4.000 


(105 codes) 


R - 


5/2 


soit 2.500 


(82 codes) 


R = 


6/3 OU 2/1 


soit 2.000 


(203 codes) 
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(suite) 
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R = 7/4 

R = 8/5 

R = 9/6 ou 3/2 

R= 10/7 

R = 11/8 

R= 12/9 ou 4/3 

R= 13/10 

R= 14/11 

R= 15/120U5/4 

R= 16/13 

R= 17/14 

R= 18/15 ou 6/5 

R = 19/16 

R = 20/17 

R = 21/18 ou 7/6 

R = 24/21 ou 8/7 



soit 1.750 
SOit 1.600 
soit 1 .500 
soit 1.428 
SOit 1.375 
SOit 1.333 
soit 1.300 
SOit 1.273 
SOit 1.250 
SOit 1.231 
SOit 1.214 
SOit 1.200 
soit 1.187 
SOit 1.176 
SOit 1.167 
SOit 1.143 



(55 codes) 
(47 codes) 
(124 codes) 
(34 codes) 
(28 codes) 
(84 codes) 
(20 codes) 
(17 codes) 
(56 codes) 
(11 codes) 
(7 codes) 
(34 codes) 
(3 codes) 
(2 codes) 
(17 codes) 
(3 codes) 
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SL^^Z 8 P ° Ur Va,6Ur ^ * N injure oo e 9 a,e a 500, . peut y avoir jusq ,a „ codes 
N«<kMM>M3dn»L. ^ m8n, "" M, " b " arel »™ 4 «''»«>«»sbM tes chac« n delon 9usu , 

ANNEXE 1 
for(x=0; x<n-m + l; X ++){ 
pos[x][0] = 0; 

for(pos[xJ[1]=p OsM [0M; pos M [1]<M-1; pos[x][1] ++ ) { 

for(pos[x][2]=pos[x][1] + l; pos[x][2]<M; pos[x]f2] ++ ) { 
(si les conditions ne sont pas satisfaites, 
continuer 
sinon, aller a ok) 

} 

} 

(arreter : impossible de Irouver un choix convenable 
pourpos[x][0...t-1]) 



ok 



} 
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ANNEXE 2 

(langage C : I'operateur ,,A " correspond au OU EXCLUSIF. 
/* Initialisation du passe de b, de longueur M 7 
for(i=M-k2; i<M; 
b[i] = 0; 

/* Calcul de N valeurs successives de la 
reponse impulsionnelle de 1/A(z) 7 
b[0] = 1; 

for(i=1; kk2; i++) 

b[i] = b[(i+M-(k2-k1)) % M] a b[(i+M-k2) % M]; 
for(i=k2; i<M; i++) 

b[i] = b[i-(k2-k1)]Ab[i-k2]; 

/* S'arranger pour qu'il n'y ait qu'un 1 

dans les k2 dernieres positions de bfiltre par A(z) 7 

weight = 0; /* tout sauf 1 7 

while(weight != 1) { 



r Decaler d'un cran 7 
for(i=1; i<M; i++) 

b[i-1] = b[ft 
b[M-1] = b[M-1-(k2-k1)] A b[M-1-k2]; 
/* Verifier 7 
weight = 0; 
for(i=M-k2; i<M; i++) { 

char sum = b[i] a b[(i+k1 ) % M] a b[(i+k2) % M]; 
if(sum) { 

shift = M - i; 

so weight++; 

} 

} 

55 
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/* Cas particulier ou M est interdit 7 
if(weight == 0) 

return(M^FORBIDDEN)- 

} 

/* decalage circulaire final a droite : 
b[ij = b[(i - shift) % MJ 7 
for(dec=0; dec < shift; dec++) { 

char temp = b[M-1]; 

for(i=M-1; i>0; i~) 
b[iJ = b[i-1J; 

b[0J a temp; 

} 



25 return(OK); 
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ANNEXE 3 

(langage C) : 

/* Phase 1 : calcul des M parites intermediaires z. 

Ces parites sont calculees en lisant les colonnes successives de la 

matrice de codage, a savoir A[*][M] A[*][N] 

Elles sont mises en tete de x a titre temporaire 7 

#define z x 
for(i=0; i<M; i++) 
= 0; 

/* Boucle sur les n-m sous-matrices a la droite de A 7 
cO = M; 
d = cO + P; 

for(k = 1 ; k <= n - m; k++) { 
offset = 0; 

for(c = c0;c<c1;c++){ 
if(x[c] !=0) 

for(i=0; kt; i++) { 

/* p devrait etre offset + pos[i]. 

On le decremente de k2 pour eviter 

le decalage 

du tableau z avant la phase 3 7 
p = offset + pos[k*t + i] -k2; 
if (p < 0) 

z[p + M] = z[p + M] A 1; 

else 

if(p < M) 

z[p] = z[p] A i; 

else 

zfp - M] = z[p - M] A 1 ; 



.109323 tAlJ_> 



11 



EP 1 093 231 A1 



10 



20 



40 



45 



} 

offset = offset + rrv 

} 

c0 = d; 
c1 =c1 +P; 

} 



/* Phase 2 : calcul des k 2 derniers symboles de parite 7 
,5 ixbO=M - 1 -k2; 

/* 1 : initialisation des k2 derniers elements de y 

temp[0...k2-1] = y[ M-1,M-2,...M-k2]7 
for(k=0; k<k2; k++) 

temp[k] = 0; 



/ 2 : recopie des k2 derniers elements de b 
reg[0...k2-1] = b[M-k2...M-1]V 
for(i=0; i<k2; i++) 

30 reg[i] = finBfi]; 

r 3 : calcul iteratif des k2 derniers symboles 7 
« for(i=0; i < M; i++) { 

/* b[i] = {l 0 0...}Ab[j.( k2 . k1)] , b[i . k2] 

avec b[i-k2]...b[i-l] = reg[0...t2-1J 

On doit verifier : 
b[-k2] + b[k1-k2] + b[0] = o 



b[-2] + b[k1-2] + b[k2-2] = 0 
b[-1] + b[k1-1] + b[k2-1] = 0 
b[0] + b[kl] + b(k2]=i 

50 b[1] + b[1+k1] + b[1+k2] = 0 

... 7 

if(i == k2) 

55 input = 1; 
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else 

input = 0; 
bi = input A reg[0] A reg[k1]; 
for(k=1;k<k2; k++) 

reg[k- 1] = reg[k]; 
reg[k2-1] = bi; 



if(bi != 0) 

is for(k=0; k<k2; k++) 
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if (z[(ixb0 - k + M) % M] != 0) 
temp[k) = temp[k] A 1 ; 

ixbO = ixbO + 1 ; 
if(ixbO = M) 
ixbO = 0; 



} 



/* 4 : Les z ont deja ete decales a gauche pour 
eviter I'ecrasement. Sinon, il faudrait faire : 

for(k=0; k<M - k2; k++) 
z[k] = z[k + k2]; 

Recopie de temp a la fin de y V 
#define y x 
for(k=0; k<k2; k++) 

y[M-1 -k] = temp[k]; 



/* Phase 3 : calcul de y[M-k2-1 , M-k2-2 o] 

y[k + k2 - k2] + y[k + k2 - k1] + y[k + k2 - 0] 4- z[k + k2] = 0 
y[k] va en x[k] 
z[k + k2] est en x[k] 
so Done : 

x[k + k2 - k2] + x[k + k2 - k1] + x[k + k2 -0] + x[k] = 0 
Soit : 

x[k + k2 - k2] = -(x[k + k2 - k1] + x[k + k2) + x[k]) 
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7 

for(k = M-k2-1;k>= 0;k~) 

yM = y[k + k2 . k1] a y[ k + k2J a Z [k]; 
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